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ABSTRACT 

We present new stellar population models of Lick absorption-line indices with variable el- 
ement abundance ratios. The models are based on our new calibrations of absorption-line 
indices with stellar parameters derived from the MILES stellar library. The key novelty com- 
pared to our previous models is that they are now available at the higher spectral resolution 
of MILES (~ 2.7 A FWHM) and flux-calibrated, hence not tied anymore to the Lick/IDS 
system. This is essential for the interpretation of galaxy spectra where calibration stars are not 
available, such as large galaxy redshift surveys or other high-redshift observations. We note 
that the MILES resolution appears to be comparable to SDSS resolution, so that our models 
can be applied to SDSS data without any corrections for instrumental spectral resolution. For 
the first time we provide random errors for the model predictions based on the uncertainties in 
the calibration functions and the underlying stellar parameter estimates. We show that random 
errors are small except at the edges of the parameter space (high/low metallicities and young 
ages < 1 Gyr) where the stellar library is under-sampled. We calibrate the base model for the 
parameters age, metallicity and a /Fe ratio with galactic globular cluster and galaxy gradient 
data. We discuss two model flavours with different input stellar evolutionary tracks from the 
Frascati and Padova groups. The new model release now includes abundance variations of the 
elements C, N, Mg, Na, Si, Ca, Ti, Cr, and Fe. The individual elements that are best accessible 
with these models and the standard set of Lick absorption features are C, N, Mg, Ca, Ti, and 
Fe. The model data is available at www.icg.port.ac. uk/^thomasd. 

Key words: stars: abundances Galaxy: abundances globular clusters: general galaxies: 
abundances galaxies: stellar content galaxies: star clusters 



1 INTRODUCTION 

The spectra of galaxies and globular clusters carry a wealth of infor- 
mation about gas and stellar population properties. Emission lines 
are used to derive gas kinematics, star formation activity and black 
hole accretion (e .g. [Ka uffmann et al. 2003; Tremont i et al.|2004| 
|Sarzi et al.|2006||Schawinski et al.|2007| >. The absorption line and 
stellar continuum component of the spectrum, instead, discloses 
stellar population properties that act as fossil record and can be 
used to derive ages, formation epochs, star formation histories, and 
element abundance s (e.g. |Trager et~a P. 2000a; Kuntschner||2000| 
Thomas et al.|2005||Nelan et al.|20051|Bernardi et al.|2006||Thomas| 



et al.|2010) . The absorption features in a spectrum are particularly 



useful for several reasons. They can easily be measured and cal- 
ibrated to a common system (Burst ein et al.| [l984; Faber et al.| 



1985 1, they are largely insensitive to dust attenuation (MacArthurl 



2005 1, and they allow to dissect 'metallicity' into individual ele- 



ment a bundances jGreggio|1997(|Tantalo et al.|19 98 ; Trag er et al.| 
|2000b[|Thomas et al.|2003) that in turn set valuable constraints on 
the chemical enrichment history (Thomas et al.l 1999 >. 



The Lick group have defined a set of 25 optical absorption- 
line i ndices jBurstein et al.|1984||Faber et al.|1985j|Gorgas et al. 



19931 |Worthey et al.||1994| |Trager et al.||1998| [Worthey & Otta 



viani |1997 1, the so-called Lick index system, that are by far the 
most commonly used in absorption-line analyses of old stellar pop- 
ulations. The bandpasses of Lick indices are relatively large with 
widths up to 50 A and two windows blue- and redward of the 
band-pass defining pseudo-continua. They have been designed for 
massive galaxies with significant line broadening through random 
motions of the stars. The advantage is that a wide bandpass in- 
creases the signal-to-noise ratios and measurements are relatively 
robust. The downside is that Lick indices are an agglomeration of a 
large number of absorption lines from various chemical elements, 
so that their use for the derivation of individual element abundances 
is nontrivial. 



|Tripicco & Be ll ( 1995 ) made a critical step forward and deter- 
mined the sensitivity of Lick absorption-line indices to individual 
element abundance variations through model atmosphere calcula- 
tions. iTTageTetal. ( 2000b ) developed a method to incorporate these 
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results in the analysis of stellar populations, which formed the bases 
for the first stellar population models of Lick absorption-line in- 
dices with variable element abundance ratios (Thomas et al. 2003 ). 
These models have subsequently been updated with new model at- 
mosphere calculations by |Korn et al.|(2~()05 see also Houdashelt et 
al 2005) for non-solar metallicities and the higher-order Balmer 
line indices ( |Thomas et al.||2004[ hereafter TMB/K models). A 
number of other element abundance sensitive models have been 
published since either for absorption-line indices following this 
semi-empirical method (Annibali et al. 2007; Schiavon 2007; Lee 
|et al.|2009| > or full spectral energy distributions using a purely the- 
oretical approach jCoelho et al. 2007 ; Wal cher et al.|2 009 1. 

Most of the models quoted above, including TMB/K, use the 
empirical calibrations of the Lick absorption-line indices with stel- 
lar parameters (the so-called fitting functions) based on the Lick 
stellar library (Worthey et al. 1994). This library has now been su- 
perseded by modern samples, most notably the Medium-resolution 
Isaac Newton Telescope library of empirical spectra ( |Sanchez^| 
|Blazquez et al.|2006| MILES). MILES features a better coverage 
in stellar parameters, higher spectral resolution, generally higher 
signal-to-noise, and a careful flux calibration. An update of the 
Maraston ( 2005 , hereafter M05) stellar population model of full 
spectral energy distributions including the MILES library is pub- 
lished in Maraston & Stromback (2010). The aim of the present 
work is to include the MILES library in the TMB/K model. To this 
end we have constructed new empirical fitting functions based on 
the MILES library {Johansson et al.|2010| hereafter JTM10). Here 
we present the final stellar population model of absorption-line in- 
dices. 

This model is the flux-calibrated version of the TMB/K model. 
It has been extended further by considering both MILES and Lick 
spectral resolution, the inclusion of additional computations based 
on the Padova stellar evolutionary tracks and a large range of indi- 
vidual element abundance variations considering the complete set 
of ten elements provided in Korn et al. (2005 ). We test this model 
on galactic globular cluster data. In a companion paper (Thomas 
|et al.|2010| hereafter Paper II) it is used to derive the element abun- 
dance ratios [C/Fe], [N/Fe], [O/Fe], [Mg/Fe], [Ca/Fe], and [Ti/Fe]. 
For the first time we provide statistical errors for each model predic- 
tion.The flux calibration in particular eases the application of these 
models to observational data where calibrating stars are not avail- 
able, ie galaxy redshift surveys such as the Sloan Digital Sky Sur- 
vey (|York et al.|2000) or oth er high-redshift observations ([Ziegler 
et al.||20051 |Bernardi et al.||20061 |Kelson et aT1|2006| |Sanchez-| 
Blazquez et al.|2009[[Thomas et al.|20T0)|Carson & Nichol|2010^ . 

The paper is organised as follows. In Sections 2 and 3 we 
present the model. The calibration of the models with galactic glob- 
ular clusters and galaxy data is shown in Section 4 as well as the 
comparison with other models in the literature. The paper con- 
cludes with Section 5. 



2 THE BASE MODEL 

The stellar population model of absorption-line indices presented 
here is an extension of the TMB/K model, which is based on the 
evolutionary population synthesis code of Maraston ( 1998, M05) 
assuming a Salpeter initial stellar mass function. We refer the 
reader to these papers for details. In the following we briefly sum- 
marise the major ingredients to the base model. 



2.1 Stellar evolutionary tracks 

Following M05 and to ensure continuity with the TMB/K model, 
the stellar evolutionary tracks from |Cassisi et al.| ( |1997| here- 
after Cassisi) are adopted, except for the highest metallicity model 
where Padova tracks (Gi rardi et al.|[20 00 hereafter Padova) are 
used. In the present work we complement our model set with fur- 
ther calculations based on Padova tracks at high metallicities with 
[Z/H] ^ —0.33 where differences are significant between the two 
model sets. Hence two separate flavours of the model based on two 
different sets of stellar tracks (Cassisi and Padova) are provided. 

Padova tracks have been chosen as comparison set, because 
these tracks are the most widely used in most other stellar pop- 
ulation synthesis models in the literature (M05, and references 
therein), in particular in other models of Lick absorption line in- 
dices (Schiavon 2007; Vazdekis et al. 2010). Cassisi and Padova 
tracks have also been chosen because they are different in some 
key aspects such as the treatment of overshooting or the mixing- 
length parameter as summarised below. In this way, the two model 
flavours presented in this paper encompass some of the uncertain- 
ties in stellar evolution theory. A major difference that impacts on 
the present model is the fact that at high metallicities, the Padova 
Red Giant Branches are cooler than those of the Cassisi tracks 
(Fig. 9 in M05). This generally leads to stronger metal absorp- 
tion indices and weaker Balmer line absorption indices (see Sec- 
tion ^. 2} , In the following we discuss some of the key ingredients 
in stellar evolutionary tracks. 



2.2 Mass loss 

Stellar mass loss cannot be predicted by stellar tracks. This comes 
from the fact that a theory relating mass-loss rates to the ba- 
sic stellar parameters does not exist. Therefore, mass loss has to 
be parametrised and its efficiency calibrated with data. In M05 
the amount of mass loss is usually parametrised by means of the 
|Reimers| {1975} empirical formula that includes a mass loss effi- 
ciency parameter r\. In order to trace the Horizontal Branch evo- 
lution properly, M05 uses the evolutionary track for the helium- 
burning phase of the mass that is obtained after mass loss is ap- 
plied to the Red Giant Branch track. Here we follow the approach 
by Maraston & Thomas ( 2000) that was to modify the mass loss 
parameter r/ as a function of metallicity such that the Balmer ab- 
sorption indices of globular clusters can be reproduced (see Sec- 
tion |4.1.1) . 



2.3 Convective overshooting 

Different from the Padova tracks, the Cassisi tracks are canonical 
stellar evolutionary tracks, i.e. the efficiency of the overshooting 
parameter is assumed to be zero (M05). The inclusion of this effect 
in stellar evolutionary tracks is still controversially discussed in the 
literature, and a detailed assessment of this issue goes beyond the 
scope of the present work. 

It is interesting to note, however, that the Padova tracks have 
a turnoff mass at given age that is larger than that of canonical 
tracks as a consequence of the different treatment of convective 
overshooting (M05). Stellar models with overshooting have more 
massive convective cores, therefore they run to higher luminosi- 
ties and live longer than classical models. Therefore, the Red Giant 
Branch phase transition occurs at 0.5 Gyr in classical models and 
at 1 Gyr in models with overshooting (M05). 
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2.4 Mixing length 

The calibration of the mixing length for the Cassisi tracks is de- 
scribed in Salaris & Cassisi ( 1996 ). The tracks with solar metallic - 
ity are computed for a = 2.25, a value that matches the Sun. This 
same value is kept in the tracks with supersolar metallicities (Bono 
|et al.|1997) . At subsolar metallicities, however, the mixing-length 
parameter is not assumed to be the same, but to vary with Z, such 
that the temperatures of the Red Giant Branch tips of Milky Way 
globular clusters are reproduced (M05). The values range from 2 
to 1.75, with a trend of decreasing mixing length with metallicity. 
In the Padova tracks, instead, the same value of the mixing-length 
parameter that is calibrated with the Sun is assumed at all metallic- 
ities. 

2.5 Non-solar a/Fe ratios 

Note that in TMB/K solar-scaled tracks were adopted, as a/Fe- 
enhanced stellar tracks at the time appeared too blue and led to 
unrealistically old ages for galaxies (Thomas & Maraston 2003]). It 
turned out that indeed the effect of a/Fe enhancement on the stel- 
lar evolutionary tracks was overestimated previously (Weiss et al. 
|2006[ >, and the recent work by Dotter et alT] {2007| > confirms that 
a/Fe enhanced stellar evolutionary tracks are not significantly dif- 
ferent from their solar-scaled counterparts. We therefore continue 
to use the solar-scaled stellar tracks employed in TMB/K. 

2.6 Stellar atmosphere 

It should be emphasised that stellar atmosphere calculations enter 
the present model only marginally. Continuum fluxes are adopted 
from the model atmospheres in order to calculate the final absorp- 
tion index. The spectral features themselves, however, that are pro- 
duced in the stars atmospheres (measured as Lick indices) are in- 
cluded the model through empirical relationships as explained in 
the following section. 



3 LICK ABSORPTION LINE INDICES 

Substantial progress in exploiting galaxy spectra has been made by 
the Lick group who defined a set of so-called absorption-line in- 
dices tracing the most prominent absorption features in the visual 
wavelength region (see Introduction). The TMB/K model has been 
developed to compute Lick absorption-line index strengths for var- 
ious chemical element abundance ratios. In the following we pro- 
vide a brief overview of the TMB/K approach with particular em- 
phasis on the novelties of the present model calculations. For more 
details on the model we refer the reader to |Thomas et al.| j2003| 
|2004l >. 



This is alternative to approaches in which absorption-line 
strengths are measured directly on the synthetic spectral energy 
distribution predicted by the stellar population model. Such mod- 
els can be based either again on empirical libraries of stellar spec- 
tra or on fully theoretical model atmosphere calculations. The lat- 
ter are known to suffer from incomplete line lists and continuum 
uncertainties, howe ver (|Korn et al.|200 5 ; Rodriguez-Merin o et al.| 
[20051 ICoelho et al.||2007| |Lee et al.||2009| |Walcher et al.||2009| >. 
As discussed in JTM10, the major strength of fitting functions, 
instead, lies in the fact that they allow for interpolation between 
well-populated regions of stellar parameter space which increases 
the accuracy of the model in stellar parameter space that is only 
sparsely sampled by empirical stellar libraries (see also Maraston 
|et al.|20 09l. Moreover, each absorption index or spectral feature is 
represented by an individual fitting function, which is optimised to 
best reproduce its behaviour in stellar parameter space. This allows 
us to compute errors in the model predictions for each individual in- 
dex (see Section |3~2} , an opportunity that is exploited in the present 
work for the first time. Finally, fitting functions are also easier to 
implement in a stellar population synthesis code, and models based 
on fitting functions are better comparable. 

In JTM10 we have recently calculated new empirical fitting 
functions for the 25 optical Lick absorption-line indices based on 
the new stellar library MILES. The MILES library consists of 985 
stars selected to produce a sample with extensive stellar parameter 
coverage. The MILES library was also chosen because it has been 
carefully flux-calibrated, making standard star-derived offsets un- 
necessary. In JTM10 we find the index measurements of the MILES 
spectra to have very high quality in terms of observational index 
errors. The new model presented here uses the JTM10 fitting func- 
tions. 



3.2 Calculation of statistical errors 

The construction of the model through fitting functions allows us 
to make a straightforward assessment of the statistical errors on 
each individual index prediction that is caused by the uncertain- 
ties in the index calibrations. We calculate errors in the model pre- 
dictions through Monte Carlo simulations. In 600 realisations per 
simulation we perturb both index measurements and stellar param- 
eters for each star using the errors given in Sanchez-Blazque z et al.| 
(2006). For each realisation we re-derive the fitting functions and 
insert those in the the stellar population code. This yields 600 index 
predictions per index and stellar population parameter age, metal- 
licity, and a/Fe ratio. We then fit a Gaussian to the distribution of 
index strengths and derive the 1-cr error. These numbers are pro- 
vided in a separate table in the model data release. Note that these 
errors do not include systematic effects such as a change of stellar 
evolutionary track, for instance. 



3.1 Empirical calibrations - fitting functions 



3.3 Spectral resolution 



Lick absorption-line strengths are computed through empirical cal- 
ibrations of index strength with stellar parameters. In this most 
widely adopted approach first absorption-line strengths are mea- 
sured on the stellar spectra of an empirical library, then analytical 
calibrations of these index strengths with stellar parameters are de- 
rived (fitting functions), which are then included in the evolutionary 
population synthesis code. In TMB/K these so-called fitting func- 
tions have been adopted from Worthey et al. ( 1994) based on the 
Lick/IDS stellar library. 



In JTM10 we have computed fitting functions for both the reso- 
lution of the MILES library and the resolution of the Lick/IDS 
library. In this work we therefore provide index predictions for 
both spectral resolutions. Note that the stellar population spectral 
resolution of MILES turns out to be somewhat worse than stated 
in |Sanchez-Blazquez etaL] (2006| l and |Vazdekis et al.| ( [2"0T0l>. We 
have run tests in which we use the GANDALF-PPXF CODE ( [Cap"^] 
Ipellari & Emsellem|2004 Sarzi et al.|2006fr and stellar population 
templates from Maraston & Stromback|i 2010) ) based on MILES 
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to derive velocity dispersions of galaxies from the Sloan Digital 
Sky Survey (SDSS) data base jYork et al.|2000} . Our results indi- 
cate that the velocity dispersions from the SDSS-MPA/JHU data 
base ([Kauffmann et al. 2003; Tremonti et al. 2004) are well re- 
produced (Thomas et al, in preparation). This indicates that stellar 
population models including the MILES library are at SDSS spec- 
tral resolution (R ~ 1800 - 2000 at 5000 A). This conclusion 
gets support from direct measurements of the spectral resolution 
on the spectral energy distributions from Mara ston & Stromback| 
( |2010} and[Vazdekis et al. (2010) m odel predictions. The latter in- 
dicate a resolution of 2.7 A FWHM |Kuntschner et al.|2006||20 10 
H. Kuntschner, private communication). 



Table 1. Element ratio bias for the various metallicities. 



Element ratio 


-2.25 


-1.35 


-0.33 


0.0 


0.35 


0.67 


[C/Fe] 


0.30 


0.30 


0.10 


0.00 


0.00 


0.00 


[N/Fe] 


0.30 


0.30 


0.10 


0.00 


0.00 


0.00 


[O/Fe] 


0.30 


0.30 


0.10 


0.00 


0.00 


0.00 


[Na/Fe] 


0.30 


0.30 


0.10 


0.00 


0.00 


0.00 


[Mg/Fe] 


0.30 


0.30 


0.10 


0.00 


0.00 


0.00 


[Si/Fe] 


0.30 


0.30 


0.10 


0.00 


0.00 


0.00 


[Ca/Fe] 


0.20 


0.20 


0.07 


0.00 


0.00 


0.00 


[Ti/Fe] 


0.15 


0.15 


0.05 


0.00 


0.00 


0.00 


[Cr/Fe] 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 



3.4 Index response functions 

Our models include different chemical mixtures and element abun- 
dance ratios. To compute optical Lick indices for these chemical 
mixtures, the impact from element ratio changes has to be as- 
sessed. This is done with the help of the so-called index response 
(2005). Extending the work of |Tripicco &| 



Korn et al. 

Korn et aL| {2005[ > calculate model atmospheres with 



functions of 
[Belli < f!995> , 

solar-scaled element ratios for various combinations of tempera- 
ture, gravity, and total metallicity. In subsequent models the abun- 
dances of the elements C, N, O, Na, Mg, Si, Ca, Ti, Fe, and Cr are 
doubled in turn, in order to determine the sensitivity of the Lick ab- 
sorption indices to element abundance variations. These response 
functions are included in the model following an extension of the 
method developed in Trager et al. (2000b). 

More specifically, Korn et al. ^2005 1 compute the model in- 
dex by splitting the basic simple stellar population model in the 
three evolutionary phases, dwarfs, turnoff stars and giants. Lick in- 
dex strengths / of the base model for each phase are computed 
separately, and modified using the fractional responses of the in- 



dex strengths SI/Iq. As discussed in Korn et al. i2005 I, these are 



applied to the flux in the absorption line rather than to the absorp- 
tion index, as the former is always a positive quantity, while indices 
can be technically negative. In this way numerical robustness is en- 
sured. 



3.5 Inclusion of element abundance variations in TMB/K 

The so-called a-elements (i.e. O, Mg, Si, Ca, Ti) and other light 
elements (i.e. C, N, Na) are combined to the 'enhanced group', 
while the iron peak elements (Fe, Cr) form the 'depressed group'. 
The ratio between those two groups is varied at fixed total metal- 
licity, and is called a/Fe. Hence the abundances of all elements in 
one group are modified in lockstep. Model predictions of Lick ab- 
sorption index strengths for the following element ratios are pro- 
vided in the TMB/K model: [a/Fe] = -0.3,0.0,0.3,0.5 dex. 
On top of this, the abundances of the individual elements carbon, 
nitrogen, and calcium are modified in separate model predictions 
providing the element ratios [a/N] = —0.5, [a/C] = —0.1, and 

[a/Ca] = -0.1, 0.2, 0.5 dex. 

The solar abundances from Grevesse et al. ( 1996 1 are adopted. 
It should be emphasised that differences between these values and 
more recent determinations can be considered small ( Asplu nd et al.| 
2009). Most importantly, however, variations of the solar abun- 
dances do not affect the models presented here. The solar abun- 
dances only serve as reference frame to which the element abun- 
dance ratios of our models are normalised. A change in solar abun- 
dance can be easily accounted for a posteriori by simple re-scaling. 



The models take into account the fact that the empirical stel- 
lar libraries used to compute model indices follow the chemical 
enrichment history of the MilkyWay, and are therefore biased to- 
wards super-solar a/Fe ratios at sub-solar metallicities. We cor- 
rect for this bias, so that the models have well-defined a/Fe ra- 
tios at all metallicities. Particular care has been taken at calibrating 
the TMB/K models with galactic globular clusters, for which ages, 
metallicities and element abundance ratios are known from inde- 
pendent sources. 



3.6 New element abundance variations 

For the new model presented here we complement and extend 
the list of element ratios provided by the TMB/K model. The 
new model is computed for the standard element ratios [a/Fe] = 
—0.3, 0.0, 0.3, 0.5 dex and additional model calculations. In each 
of these, only one of the elements C, N, Na, Mg, Si, Ca, and Ti is 
increased by 0.3 dex relative to the a-element abundance. Hence, 
this yields seven different additional model tables with the element 
ratios [X/a] = 0.3 dex for X =C, N, Na, Mg, Si, Ca, Ti. These 
can be used to assess the effect of abundance variations of these el- 
ements. We do not include O, because this element dominates total 
metallicity, and Fe, because this element is already considered with 
varying a/Fe ratio. 

We emphasise again that metallicity is kept fixed in these cal- 
culations. Note, however, that these elements do not contribute sig- 
nificantly to the overall mass budget when treated individually. 
Therefore, one can safely consider those enhancements as minor 
perturbations. Since this does not apply to O, we do not calculate 
individual enhancements for this element. As discussed in Paper II, 
there is strong evidence now that the heavier of the a elements like 
Ca and Ti are less enhanced than O or Mg in metal-poor Milky Way 
stars ( |Feltzing et al.|2009||Bensby et al.|2010) . Different from the 
TMB/K model we consider these elements separately in the new 
model, therefore this differential bias in the library stars has to be 
taken into account (see Section [33} , We assume a bias of 0.3 dex 
for the light elements, which we reduce to 0.2 — 0.15 dex for those 
elements with higher atomic number. The detailed bias adopted is 
summarised in Table [T] Note that the element ratio [O/Fe] is in- 
cluded in the bias table, even though we do not consider individ- 
ual O/Fe ratios. The reason is that O abundance is included in the 
'enhanced group' . Like in TMB/K we do not assume any bias at 
super-solar metallicity as suggested by Pro ctor et al.| ([2"004), as the 
various elements still exhibit conflicting trends (see discussion in 
|Thomas et al.|2005| >. 
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3.7 Summary of new features 

The new features of this model are. 1) Based on a flux-calibrated 
version of the Lick/IDS index system. This is achieved through 
the new flux calibrated fitting functions of JTM10 based on the 
MILES stellar library. 2) Error estimates for each index as function 
of stellar population parameters. 3) Model predictions provided for 
MILES and Lick spectral resolutions. 4) Two different stellar evo- 
lutionary tracks (Cassisi and Padova). 5) Additional sets of models. 
In each of these, only one of the elements C, N, Na, Mg, Si, Ca, 
and Ti is increased by 0.3 dex relative to the a-element abundance. 
A differential element ratio bias in the stellar libraries at low metal- 
licities is considered to account of the fact that heavier a elements 
tend to be less enhanced in metal-poor halo stars. 



4 MODEL TESTING 

Following our strategy for the TMB/K models, we compare the 
model predictions with observational data of galactic globular clus- 
ters, as the latter are the closest analogues of simple stellar popu- 
lations in the real universe (Maraston et al. 2003}. Key is that in- 
dependent estimates of ages, metallicities, and element abundance 
ratios are available for the globular clusters of the Milky Way from 
deep photometry and stellar high-resolution spectroscopy. As dis- 
cussed in Thomas et al. (2003 . see references therein) galactic glob- 
ular clusters are known to be old and a/Fe enhanced similar to halo 
field stars (for more recent compilations see |De A ngeli et al .|2005| 
|Pritzl et al.|2005l|Mendel et al.|2007> . 

The globular cluster samples considered here are from |Puzia| 
|et al.| l |2002"} and |Schiavon et al.| ( [2005} . We do not use the indices 
tabulated in |Puzia ^T al. (2002) directly, because these measure- 
ments have been calibrated onto the Lick/IDS system by correcting 
for Lick offsets. Schiavon et al.|p005| > do not provide line index 
measurements. Hence we measure line strengths for all 25 Lick 
absorption-line indices directly on the globular cluster spectra us- 
ing the definitions by Trager et al. (1998 ) and Worthey & Ottaviani| 
( |1997) . Both globular cluster sample have been flux calibrated, so 
that no further offsets need to be applied for the comparison with 
the models presented here. We have smoothed the spectra to Lick 
spectral resolution before the index measurement in order to allow 
for a straight compar ison with the TMB/K model. Note that the 



Table 2. Reduced x 2 values for model-data comparison in Figs. [T] — [4] and 
median relative errors in the globular cluster sample. 



Schiavon et al. 



1 2005 1 spectra are corrupt around 4546 and 5050 A, 
so that the indices Fe4531 and Fe5015 cannot be measured (S05). 
For more details we refer the reader to Paper II. A table with the 
resulting indices is provided in the online model data release. 

The aim of this section is to test the new models with globu- 
lar cluster data. As index definitions have not been modified in the 
present work, we refer the reader to the literature for a comprehen- 
sive discussion of Lick index behaviour with element abundance 
ratio (e.g.|Tripicco & Bell|1995)|Trager et al.|2000b[|Thomas et al.| 
|2003"l|2004||Kornetal.|2005| >. 



4.1 Comparison with globular cluster data 



Figs. [T| [2] [3] [4] are are-make of Fig. 2 from |Thomas et aL] ( |2003| >. 
Following Thomas et al. (2004) we plot predictions for the 25 
absorption-line indices as functions of the index [MgFe]' in order 
to visualise better the a/Fe sensitivity of each individual index. 
This index has been defined in Thomas et al. (2 003) as 
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in order to eliminate the residual a/Fe dependence of the index 
[MgFe] introduced by Gonzalez ( 1993 ). The plots show three mod- 
els at Lick spectral resolution with an age of 13 Gyr, the metal- 
licities [Z/H] = -2.25, -1.35, -0.33, 0.0, 0.35, 0.67 dex, 
and [a/Fe] = 0.0, 0.3, 0.5 dex as solid magenta, yellow, and cyan 
lines, respectively. Hence each line is a model at fixed age and a/Fe 
ratio with total metallicity increasing from left to right. The grey 
shaded area along the model with [a/Fe] = 0.3 dex (yellow line) 
indicates the 1-cr error of the model prediction. 

Dotted lines are the original TMB/K model for comparison. 
Note that this model is inherent to the Lick/IDS system. Discrepan- 
cies between the TMB/K and the new model will be caused mostly 
by differences in the fitting functions adopted (Worthey et al vs 
JTM10), and the flux-Lick/IDS calibration. 

Data of galactic globular clusters from Puzia et al. (2002} and 
Schia voneTaL] |2005) are filled squares and open squares, respec- 
tively. The typical errors in the globular cluster index measurements 
(see Paper II) are given by the error symbol at the bottom of each 
panel. Finally, the small black dots are early-type galaxies from the 
MOSES catalogue (Morphologically Selected Early-type galaxies 
in SDSS; [Schawinski et al.|2007l|Thomas et al.|2010} drawn from 
the SDSS data base. Only measurements on high signal-to-noise 
spectra with S/N > 40 are shown, primarily in order to provide an 
indication for the locus of massive galaxies in this diagram. 

We split the presentation in four different figures: the five 
Balmer line indices in Fig. [T](H<5a, H5f, H7A, H7F, H/3), the 
indices with a positive response to a/Fe enhancement in Fig. [5] 
(CNi, CN 2 , Ca4227, G4300, Mg l5 Mg 2 , Mgfc), the indices 
with negative responses to a/Fe enhancement in Fig.[3](Fe4383, 
Fe4531, Fe5015, Fe5270, Fe5335, Fe5406, Fe5709, Fe5782) , 
and the remaining indices that are insensitive to a/Fe in Fig. |4] 



6 Thomas, Maraston, Johansson 




Figure 1. Calibration of the Balmer line indices. Three models at Lick spectral resolution with an age of 13 Gyr, the metallicities [Z/H] = 
— 2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] = 0.0,0.3,0.5 dex are shown as solid magenta, yellow, and cyan lines, respectively. Each 
line is a model at fixed age and a/Fe ratio with total metallicity increasing from left to right. The grey shaded area along the model with [a/Pe] = 0.3 dex 
(yellow line) indicates the l-cr error of the model prediction. Dotted lines are the original TMB/K model that is inherent to the Lick/IDS system. Galactic 
globular clusters from Puzia et al. ( 2002 1 and Schiavon et al. ( 2005 1 are filled and open squares, respectively. The typical errors in the globular cluster index 
measurements are given the error symbol at the bottom of each panel. The small black dots are early-type galaxies from the MOSES catalogue (Morpholog- 
ically Selected Early-type galaxies in SDSS Schawinski et al. 2007 Thomas et al. 2010) drawn from the SDSS (Sloan Digital Sky Survey) data base {York| 
|et al.|20do") including only high signal-to-noise spectra with S/N > 40. 



(Ca4455, C 2 4668, NaD, TiOi, Ti0 2 ). Following [Thomas et al.| 
( 2003 2004 1 we consider a model well calibrated if the model track 
with [a/Fe] = 0.3 matches the globular cluster data. 

To quantify this comparison, we calculate the median reduced 
X 2 values for this comparison. We perform two tests, one including 
all globular clusters, and one focusing on high metallicity clusters 
with [Z/U] > -0.8 dex. The results are listed in Table [5] Note 
that these numbers should be interpreted with care. Their aim is to 
aid the identification of well calibrated and badly calibrated indices 
in the following discussion. Low \ 2 can be produced by large er- 
rors, and may not necessarily indicate a good match. Furthermore, 
this exercise does not include the variation of other element ratios, 
that we know are crucial for some indices, in particular those that 
are sensitive to individual elements whose abundances may show 
anomalies in galactic globular clusters. For instance, the indices 
CNi and CN 2 are very sensitive to N abundance, and it is known 
that nitrogen abundance is additionally enhanced in galactic glob- 



ular clusters (Thomas et al. 2003 see discussion in Paper II and 
references therein). A proper fit to the data including the variation 
of various other element abundance ratios is presented in Paper II. 

Therefore, we deem an index badly calibrated if the model 
prediction for [a/Fe] = 0.3 dex deviates significantly from the 
globular cluster data (chi 2 ^ 1.5 in either Columns 2 or 3 in Ta- 
ble |2j, if the median observational relative error in the globular 
cluster sample exceeds 0.5 (Column 4 in Table |2j, and if the index 
is not sensitive to other element ratios beyond the a/Fe ratio (Fig. 1 
in Paper II). The result of this selection process is summarised in 
Column 5 of Table tab:chi2. 

Finally, it is interesting to note that model errors are generally 
very small and well below the observational errors around solar 
metallicity. Errors rise considerably and become comparable to the 
typical observational error toward the highest and lowest metallic- 
ities. This behaviour ought to be expected and is a direct conse- 
quence of the empirical stellar library, in which the stellar param- 
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Figure 2. Calibration of the indices that correlate with a/Fe. Three models at Lick spectral resolution with an age of 13 Gyr, the metallicities [Z/H] = 
—2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] =0.0, 0.3, 0.5 dex are shown as solid magenta, yellow, and cyan lines, respectively. Each line is a 
model at fixed age and a/Fe ratio with total metallicity increasing from left to right. The grey shaded area along the model with [a/Fe] =0.3 dex (yellow line) 
indicates the 1-cr error of the model prediction. Dotted lines are the original TMB/K model that is inherent to the Lick/IDS system. Galactic globular clusters 
from |Puzia et al. 1^2002) and |Schiavon et aL](2005) are filled and open squares, respectively. The typical errors in the globular cluster index measurements 
are given the error symbol at the bottom of each panel. The small black dots are early-type galaxies from the MOSES catalogue (Morphologically Selected 
Early-type galaxies in SDSS |Schawinski et al.|2007||Thomas et al.|2010) drawn from the SDSS (Sloan Digital Sky Survey) data base ^York et al.|2000) 
including only high signal-to-noise spectra with S/N > 40. 



eter space is inevitably sampled worse at the ends of the distribu- 
tion. This is well consistent with a recent conclusion of IVazdekisl 
|et al.| ( [2bl0| l who show that the stellar population model quality 
decreases toward the edges of parameter space. 



4.1.1 Balmer line indices 

Fig. [T] confirms our previous results that the Balmer line indices 
are generally quite sensitive to a/Fe ratio variations. Absorption- 
line strengths of the higher-order Balmer lines H<5a, H5p, and H7A 
increase considerably with increasing a/Fe at super-solar metal- 
licities. This is critical for the derivation of galaxy ages (jThomas 
|et al.|20 04 , Thomas & Davies 2006). As discussed in Tho mas et al. 
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Figure 3. Calibration of the indices that anti-correlate with a/Fe. Three models at Lick spectral resolution with an age of 13 Gyr, the metallicities [Z/H] = 
—2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] = 0.0, 0.3, 0.5 dex are shown as solid magenta, yellow, and cyan lines, respectively. Each line is 
a model at fixed age and a/Fe ratio with total metallicity increasing from left to right. The grey shaded area along the model with [q/Fc] = 0.3 dex (yellow 
line) indicates the \-a error of the model prediction. Dotted lines are the original TMB/K model that is inherent to the Lick/IDS system. Galactic globular 
clusters from Puzia et al. 1 2002 1 and Schi avon et al.H2005) are filled and open squares, respectively. Note that we could not make reliable measurements of 
Fe4531 and Fe5015 index strengths for the |Schiavon et al.| (2005') sample. The typical errors in the globular cluster index measurements are given the error 
symbol at the bottom of each panel. The small black dots are early-type galaxies from the MOSES catalogue (Morphologically Selected Early-type galaxies 
in SDSS [Schawinski et al.|2007||Thomas et al.|20"T0) drawn from the SDSS (Sloan Digital Sky Survey) data base {York et al.(2000) including only high 
signal-to-noise spectra with S/N > 40. 



(2004 ) and Korn et al. (2005 ) this is caused by prominent Fe ab- 
sorption in the pseudo-continuum windows of the index definition. 
The indices H7F and H/3, instead, are only very little sensitive to 
a/Fe. 

The globular cluster data seems generally well reproduced by 
the model for all B aimer line indices. In particular the trend of in- 
creasing Balmer line strength with decreasing metallicity is well 



matched. In more detail, however, it can be seen that globular clus- 
ters at higher metallicities are below the models for the indices H<5p 
and, most prominently, H/5. The \ 2 values quoted in Table[2]con- 
firm this conclusion. Other recent models in the literature show 
the same patte rn (|Schiavon||2007| |Lee et al.||2009| [Cervantes &| 
|Vazdekis|2009l|Vazdekis et al.|2010[|Poole et al.|2010> , which has 
therefore been called 'H/3 anomaly' by Poole et al. (2010). Fig. |TJ 
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Figure 4. Calibration of the indices that are insensitive to a/Fe. Three models at Lick spectral resolution with an age of 13 Gyr, the metallicities [Z/H] = 
— 2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] = 0.0, 0.3, 0.5 dex are shown as solid magenta, yellow, and cyan lines, respectively. Each line is a 
model at fixed age and o/Fe ratio with total metallicity increasing from left to right. The grey shaded area along the model with [a /Fe] =0.3 dex (yellow line) 
indicates the 1-cr error of the model prediction. Dotted lines are the original TMB/K model that is inherent to the Lick/IDS system. Galactic globular clusters 
from |Puzia e7aE1{2002) and |Schiavon et al.|(2005) are filled and open squares, respectively. The typical errors in the globular cluster index measurements 
are given the error symbol at the bottom of each panel. The small black dots are early-type galaxies from the MOSES catalogue (Morphologically Selected 
Early-type galaxies in SDSSjSchawinski et al. 2007 Thomas et al. 2010) drawn from the SDSS (Sloan Digital Sky Survey) data base (York et al. 2000) 
including only high signal-to-noise spectra with S/N > 40. 



shows that this anomaly may extend to the other Balmer lines, at 
least H<5f, while the other higher-order Balmer line indices seem 
not to be affected. The location of the index measurements for 
SDSS galaxies is well covered by the models. Younger ages would 
need to be considered, of course, to match objects with stronger 
Balmer line indices. 

The TMB/K model is shown by the dotted lines. Offsets be- 
tween Lick and flux-calibrated systems are clearly small, as the pre- 
dictions from the new model are almost identical to TMB/K. The 
significantly lower Balmer line strengths at low metallicities come 
partly from lower index values in the fitting functions for hot gi- 
ant stars (JTM10) and partly from a re-adjustment of the mass loss 
parameter r\ such that the globular cluster data are matched. The 
efficiency of mass loss along the Red Giant Branch evolution is ex- 
pected to decrease with deceasing metallicity (e.g. |Renzini & Volt] 
1 1 98 1 [ >. To minimise the problem of the Balmer anomaly, we had to 
set r) = at the lowest two metallicities. Note that the Horizontal 



Branch morphology is still not completely red, because the metal- 
licity is so low that the Horizontal Branch track has naturally a high 
enough effective temperature such that a non-negligible amount of 
fuel is spent bluewards of the RR Lyrae strip even without mass 
loss (Fig. 1 1 in M05). 



4.1.2 Indices that correlate with a/Fe 

Fig. [2] presents the calibration of the indices whose strengths in- 
crease with an increase of the a/Fe ratio, i.e. CNi, CN2, Ca4227, 
G4300, Mg x , Mg 2 , and Mg6. Note that most of this correlation 
with a/Fe is actually caused by an anti-correlation with Fe abun- 
dance (Thomas et al. 2003 1. The index that is most sensitive to the 
a /Fe ratio is Mg b. Globular cluster data are very well matched 
for G4300, Mg 2 , and Mgb. The models are well off, instead, for 
the indices CNi, CN 2 , Ca4227 as confirmed by the high \ 2 val- 
ues in Table [2] Mgj^ appears to be slightly below the model, but 
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Figure 5. Calibration of Balmer line indices, equivalent to Fig. ^ Three models at Lick spectral resolution with an age of 13 , the metallicities [Z/H] = 
—2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] = 0.3 dex are shown. The yellow lines are the same model as in Fig. phased on the Cassisi stellar 
evolutionary tracks. The grey shaded area along the model indicates the 1-cr error of the model prediction. The magenta lines are the model based on Padova 
tracks at [Z/H] ^ —0.33 dex. The cyan lines are the m odel by|Schiavon| 2007}. Each line is a model at fixed age and a/Fe ratio with total metallicity 
increasing from left to right. Galactic globular clusters from Puzia et al. 1 2002 and Schiavon et al. ( 2005! are filled and open squares, respectively. The typical 
errors in the globular cluster index measurements are shown as the error symbol at the bottom of each panel. The small black dots are early-type galaxies from 
the MOSES catalogue (Morphologically Selected Early-type galaxies in SDSS Schawinski et al. 2007, Thomas et al. 2010 1 drawn from the SDSS (Sloan 
Digital Sky Survey) data base jYork et al.|2000} including only high signal-to-noise spectra with S/N > 40. 



the observational error is large enough to keep the \ 2 l° w - This 
is caused by the variation of further chemical elements these in- 
dices are sensitive to. As already discussed in Thomas et al. (2003 1, 
the stronger CN indices in the globular cluster data indicate an en- 
hancement of nitrogen, while the weaker indices Ca4227 and Mgj^ 
can be explained by a depression of calcium and carbon. In Paper II 
we present the full analysis of these element abundance variations 
and show that the globular cluster data can be recovered very well. 

It can further be seen that SDSS galaxies are on line with with 
globular clusters as far as the enhancement of the a/Fe ratio is 
concerned, but may well have different chemical mixtures regard- 
ing other elements like carbon, nitrogen, or calcium. A full chemi- 
cal analysis of SDSS galaxies is the subject of a companion paper 
(Johansson et al, in preparation). 

Finally, the TMB/K model (dotted lines) is quite consistent 
with our new, flux-calibrated model for the indices Ca4227, G4300, 
and Mg b. Deviations are significant, instead, for CNi, CN2, Mg ls 



and Mg 2 , even though the fitting functions used in TMB/K and in 
this work are consistent for these indices and do not explain this 
offset (JTM10). The most likely cause for deviations is the fact that 
these indices have their pseudo-continuum windows placed rela- 
tively far from the actual index bandpass, so that variations of the 
shape of the stellar continuum play a larger role. 



4.1.3 Indices that anti-correlate with a/Fe 

Fig. [3] presents the calibration of the indices whose strengths de- 
crease with an increase of the a/Fe ratio, i.e. all Fe indices Fe4383, 
Fe4531, Fe5015, Fe5270, Fe5335, Fe5406, Fe5709, and Fe5782. 
Clearly, the anti-correlation with a/Fe comes from the response 
to the depression of Fe abundance. The indices most sensitive to 
a/Fe are Fe4383, Fe5335, and Fe5406. The globular cluster data 
are reproduced well for indices Fe4383, Fe4531, Fe 5270, Fe5335, 
and Fe5406 as also indicated by the low x 2 m Table u\ The mod- 



Stellar population models of Lick absorption-line indices 1 1 



TT — I — I | I I I I — | I I I I | — I I I I | I I I I | I I I I I I I — | I I I I | I I I I | I I I I — | I I I I — | I I 




Figure 6. Comparison of the indices that correlate with ct/Fe, equivalent to Fig. [2] Three models at Lick spectral resolution with an age of 13 , the metallicities 
[Z/H] = —2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] = 0.3 dex are shown. The yellow lines are the same model as in Fig.[2]based on the 
Cassisi stellar evolutionary tracks. The grey shaded area along the model indicates the 1-cr error of the model prediction. The magenta lines are the model 
based on Padova tracks at [Z/H] ^ —0.33 dex. The cyan lines are the model by|SchiavonM2007k. Each line is a model at fixed age and a/Fe ratio with total 
metallicity increasing from left to right. Galactic globular clusters from |Puzia et al.H 2002 1 and Schiavon et al. 1 2005 1 are filled and open squares, respectively. 
The typical errors in the globular cluster index measurements are shown as the error symbol at the bottom of each panel. The small black dots are early-type 
galaxies from the MOSES catalogue (Morphologically Selected Early-type galaxies in SDSS[Schawinski et al. 2007, Thomas et al. 2010) drawn from the 
SDSS (Sloan Digital Sky Survey) data base I York et al. 2000) including only high signal-to-noise spectra with S/N > 40. 



els are clearly off the data for Fe5015 and Fe5782, which is further 
supported by the high \ 2 values. Fe5709 has a low \ 2 m spite of a 
bad match, which is caused by the relatively large measurement er- 
rors on this index (see Table|2|. For the badly calibrated indices, is 
not clear at this point whether the problem lies in the model or the 
globular cluster data, but it is certainly advisable not to use these 
two indices in an analysis of absorption line spectra. 

Most model predictions are consistent with the TMB/K model. 



The new model predicts lower index strengths at high metallicities 
for the indices Fe4383, Fe4531, and Fe5015, which is almost en- 
tirely caused by the offset between the Lick and the flux-calibratd 
systems. 
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Figure 7. Comparison of the indices that anti-correlate with a/Fe, equivalent to Fig. [3] Three models at Lick spectral resolution with an age of 13 , the 
metallicities [Z/H] = —2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] = 0.3 dex are shown. The yellow lines are the same model as in Fig. [5] 
based on the Cassisi stellar evolutionary tracks. The grey shaded area along the model indicates the 1-cr error of the model prediction. The magenta lines are 
the model based on Padova tracks at [Z/H] ^ —0.33 dex. The cyan lines are the model by Schiavon|12007k. Each line is a model at fixed age and a/Fe ratio 
with total metallicity increasing from left to right. Galactic globular clusters from |Puzia et al.| j2002l and Schiavon et al. (2005 1 are filled and open squares, 
respectively. Note that we could not make reliable measurements of Fe4531 and Fe5015 index strengths for the [Schiavon et al. (20051 sample. The typical 
errors in the globular cluster index measurements are shown as the error symbol at the bottom of each panel. The small black dots are early-type galaxies from 
the MOSES catalogue (Morphologically Selected Early-type galaxies in SDSS Schawinski et al. 2007, Thomas et al. 20101 drawn from the SDSS (Sloan 
Digital Sky Survey) data base jYork et al.|2000} including only high signal-to-noise spectra with S/N > 40. 



4.1.4 Indices insensitive to a/Fe 

Finally, the indices that are not sensitive to a/Fe ratio variations 
are shown in Fig.|4]. These are Ca4455, C 2 4668, NaD, TiOi, and 
Ti02- The three reddest indices are clearly not well calibrated as 
has already been concluded in Thomas et al. (2003 1, and we can- 
not recommend these features for the analysis of absorption line 



spectra. Again, the relatively large observational errors in the TiO 
index measurements lead to low x 2 values despite the clear offset 
between models and observations (see Table|2]l. NaD is notoriously 
difficult because it is affected by interstellar absorption, so that the 
observational measurement becomes meaningless. The same dif- 
ficulty seems to apply to TiOi, in which case the interstellar Na 
absorption affects the blue pseudo-continuum window and leads to 
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Figure 8. Comparison of the indices that are insensitive to a/Fe, equivalent to Fig. [4] Three models at Lick spectral resolution with an age of 13 , the 
metallicities [Z/H] = —2.25, —1.35, —0.33, 0.0, 0.35, 0.67 dex, and [a/Fe] = 0.3 dex are shown. The yellow lines are the same model as in Fig. [4] 
based on the Cassisi stellar evolutionary tracks. The grey shaded area along the model indicates the 1-cr error of the model prediction. The magenta lines are 
the model based on Padova tracks at [Z/H] ^ —0.33 dex. The cyan lines are the model by|Schiavon |2007}. Each line is a model at fixed age and a/Fe ratio 
with total metallicity increasing from left to right. Galactic globular clusters from Puzia et al. (2002 and Schiavon et al. (2005! are filled and open squares, 
respectively. The typical errors in the globular cluster index measurements are shown as the error symbol at the bottom of each panel. The small black dots are 
early-type galaxies from the MOSES catalogue (Morphologically Selected Early-type galaxies in SDSS Schawinski et al. 2007 Thomas et al. 2010) drawn 
from the SDSS (Sloan Digital Sky Survey) data base {York et al.|2000} including only high signal-to-noise spectra with S/N > 40. 



weaker index measurements. Note also that the scatter in the lSchi-l 
avon (2007 ) cluster measurements is unusually large for all three 
indices, which further hampers a meaningful model calibration. It 
is less clear, though, why the model fails for TiC>2. 



Ca4455 is reasonably well calibrated. Most notably, the cali- 
bration has improved significantly. The TMB/K model seemed to 
fail in reproducing observations, most likely because of a wrong 
Lick calibration of the ] 



Puzia et al. 1 2002 1 data as discussed already 



in Maraston et al. (20031 and Thomas et al. (2003). This problem 
has been solved now with our new flux-calibrated models. Still, 
we exclude Ca4455 from the recommended set of indices because 
of its relatively large measurement error (see Table [2}. C24668 is 
offset as supported by the enhanced x 2 in Table [2] This particular 
index is sensitive to C abundance, though, so the mismatch might 
well be an abundance ratio effect. 



4.2 Models with Padova tracks and literature comparison 

In this section we compare with our second set of models based 
on the Padova stellar evolutionary tracks at high metallicities 
([Z/R] ^ -0.33 dex) and with the models of|Schiavon|(|2007^. As 



far as we are aware, the latter is the only other model in the current 
literature that includes variable element ratios and is flux-calibrated 
hence not tied to the Lick/IDS system. The figures showing this 
comparison equivalent to Figs.[TJ[2]|3] and|4]are Figs.|5]|6]|7] and[8] 
The plots show again three models at Lick spectral resolution with 
an age of 13 Gyr, all six metallicities, and a/Fe = 0.3 dex. The 
yellow lines are the same model as in Fig. [TJ based on the Cassisi 
stellar evolutionary tracks. The magenta lines are the model based 
on Padova tracks at high metallicities, and the cyan lines are the 
model of |Schiavon| l |200 ! 7) , 

There is good agreement between our model and the model 
by |Schiavon| ( |2007l l for the indices Mg 2 , Mg6, Fe5015, Fe5270, 
and Fe5335. The higher-order Balmer line strengths tend to be 
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Figure 9. Remake of Fig. 3 in |Kuntschner et al.|{2010) . H/3 is plotted as a function of [MgFe50]' {Kuntschner et al.|2010| . The SAURON central galaxy 
data are the circles with the tails indicating the index gradients to larger radii. Stellar population models are over-plotted for various ages and metallicities at 
fixed solar a/Fc. Models based on Cassisi and Padova tracks are shown in the left-hand and right-hand panels, respectively. The age of the oldest model track 
is 13 Gyr. Different from Fig. 3 in Kuntschner et al. (20101 this plot compares observational and model data in flux-calibrated space independently of Lick 
offsets. 



somewhat higher in the Schiavon ( 2007 I model at solar and super- 
solar metallicities. The largest discrepancies are found for CNi, 
CN 2 , Ca4227, G4300, Fe4383, and C 2 4668. Note that in all these 
cases the discrepancies are far larger than the differences between 
our models with different evolutionary tracks (yellow and magenta 
lines). 

Our model produces significantly stronger CNi and CN 2 at 
all metallicities, which is better in line with the globular cluster 
data. Likewise, the slightly higher Fe4383 index prediction in our 
model gets support from the observations. The G4300 index also is 
lower in the Schiavon ( 2007 I models, but fits the data better than 
our model. Ca4227 is significantly stronger in |Schiavon| (2007). 
The comparison with the globular cluster data is less decisive in this 
case, however, as both model deviations can be corrected through 
a modification of Ca abundance. The Schiavon] ( |2007| > models 
will lead to considerably lower Ca abundances. Finally, C 2 4668 is 
stronger in our models only at high metallicities, which inhibits the 
direct comparison with globular cluster data. The Schiavon ( 2007 I 
are expected to yield larger carbon abundances for galaxies. 

Our model based on Padova tracks (magenta lines) is consis- 
tent with the base model for the majority of the indices. The cases 
for which the discrepancy significantly exceeds the model error 
(grey shaded area) are H/3, CN X , CN 2 , and C 2 4668. Small de- 
viations are found for Ca4227, G4300, Ca4455, Fe5015, Fe5709. 
In all cases, equally for the Balmer and the metal lines, the Padova 
based models produce lower index strengths. This effects kicks in 
at super-solar metallicities, however, where no globular cluster data 
is available for calibration. 



4.3 Comparison with galaxy data 

We therefore turn to consider galaxy data for a more detailed com- 
parison between our models sets with Cassisi and Padova tracks. 
We use the recently published index gradient data from SAURON 
( |Kuntschner et al . 2006 2010). This data provides high signal-to- 
noise measurements of the absorption-line indices H/3, Fe5015, and 
Mg b. Most critically for the comparison with our new models is 



that this data has been carefully flux calibrated. Fig. [9]is a remake 
of Fig. 3 in Kuntschner et al. (2010) plotting H/3 as a function of 



[MgFe50]'. This latter index has been defined in Kuntschner et al 



( 2010 1 in the style of [MgFe]' combining Mg b with Fe5015 rather 



than Fe5270 and Fe5335 because of the restricted wavelength range 
sampled in SAURON observations. 



[MgFe50]' 



0.69 x Mg6 + Fe5015 



(2) 



The stellar population model grid is shown for various ages 
and metallicities (at solar a/Fe) on top of the galaxy data. Mod- 
els based on Cassisi and Padova tracks are shown in the left-hand 
and right-hand panels, respectively. The age of the oldest model 
track is 13 Gyr. Different from Fig. 3 in |Kuntschner et aT](2010| ) 
this plot compares observational and model data in flux-calibrated 
space independently of Lick offsets. We confirm the conclusion by 
|Kuntschner et al.| ( |2010| l that observations at all galaxy radii are rea- 
sonably well covered by the models. While globular cluster data do 
not allow for the calibration at these high metallicities, the galaxy 
data clearly suggests that the model based on Padova tracks per- 
forms better. Note that the oldest model track plotted is 13 Gyr, 
hence well within the age of the universe (Komatsu et al. 20T0). 

We conclude that galaxy data do not confirm the existence of 
an H/3 anomaly indicated by globular cluster data ( |Poole et al.| 
2010). This indicates that the cause for the Balmer anomaly in 
globular clusters is more likely to be found in the observational 
measurement. The origin of the problem could be a failure of the 
observations to sample all evolutionary phases properly (Maraston 
|et al.| |2003 l. Fig. [9] suggests that here is no problem with galaxy 
age-dating through the H/3 index, instead. 



5 CONCLUSIONS 

We present new stellar population models of absorption-line in- 
dices with variable element abundance ratios. The model is an ex- 
tension of the TMB/K model, which is based on the evolutionary 
stellar population synthesis code of Maraston (1998, 2005J. The 
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new models are based on our new calibrations of absorption-line 
indices with stellar parameters derived from the MILES stellar li- 
brary (Johansson et al. 20101. The MILES library consists of 985 
stars selected to produce a sample with extensive stellar parame- 
ter coverage. The MILES library was also chosen because it has 
been carefully flux-calibrated, making standard star-derived offsets 
unnecessary. 

The key novelty compared to our previous models is that they 
are flux-calibrated, hence not tied anymore to the Lick/IDS system. 
This is essential for the interpretation of galaxy spectra where cali- 
bration stars are not available, such as large galaxy redshift surveys 
like SDSS or other high-redshift observat ions {Ziegler et al.|2005| 
Bernardi et al. 2006; Kelson et al. 2006, Sanchez-Blazquez et al. 
2009| |Thomas et al.|2010||Carson & Nichol|2010) . A further new 
feature is that model predictions are provided for both the original 
Lick and the higher MILES spectral resolutions. Note that the latter 
turns out to be comparable to the SDSS resolution, so that our new 
high-resolution models can be applied to SDSS data without any 
corrections for instrumental spectral resolution. 

The construction of the model through fitting functions allows 
us to make a straightforward assessment of the statistical errors on 
each individual index prediction. Hence as a further novelty we cal- 
culate errors in the model predictions through Monte Carlo simu- 
lations, which are provided in a separate table in the model data 
release. It turns out that the model errors are generally very small 
and well below the observational errors around solar metallicity. 
Errors rise considerably and become comparable to or may exceed 
the typical observational error toward the highest and lowest metal- 
licities. This behaviour ought to be expected and is a direct conse- 
quence of the empirical stellar library, in which stellar parameter 
space is inevitably sampled worse at the ends of the distribution. 

The data release now provides models with two different stel- 
lar evolutionary tracks by Cassisi (used in TMB/K) and Padova at 
high metallicities. The model based on Padova tracks is consistent 
with the base model for the majority of indices. The cases of indices 
for which the discrepancy exceeds the model error significantly 
are H/3, CNi, CN2, and C24668. Small deviations are found for 
Ca4227, G4300, Ca4455, Fe5015, and Fe5709. In all cases, equally 
for the Balmer and the metal lines, the Padova based models pro- 
duce lower index strengths. 

Finally, as the last novelty of this model, we release additional 
models tables with enhancement of each of the elements C, N, Na, 
Mg, Si, Ca, Ti, and Cr separately by 0.3 dex. 

We calibrate the base model for the parameters age, total 
metallicity and a/Fe ratio with galactic globular cluster and galaxy 
gradient data. Key is that independent estimates of ages, metallic- 
ities, and element abundance ratios are available for the globular 
clusters of the Milky Way from deep photometry and stellar high- 
resolution spectroscopy. The globular cluster samples considered 
here are from P uzia et al.| |2002l and Sc hiavon et al.| ( [2005| ). For 
both samples we measure line strengths of all 25 Lick absorption- 
line indices directly on the globular cluster spectra (Paper II). Both 
globular cluster samples have been flux calibrated, so that no fur- 
ther offsets need to be applied for the comparison with the models 
presented here. 

The globular cluster data is well reproduced by the model for 
all Balmer line indices. In particular the trend of increasing Balmer 
line strength with decreasing metallicity is well matched. We con- 
firm previous findings of a mild H/3 anomaly with models gener- 
ally predicting too strong H/3 indices at intermediate metallicities. 
We show that a similar problem, even though slightly less severe, 
exists for H<5f . The problem is not replicated, however, in the com- 



parison with galaxy data at similar metallicities from Kun tschner] 
|et al.| < [2010| ). A good match with globular cluster data is seen for 
the a/Fe sensitive, metallic indices G4300, Mg 2 , and Mgfo. The 
models are well off, instead, for the indices CNi, CN2, Ca4227, 
C24668, and Mgj, which is caused by the variation of further 
chemical elements these indices are sensitive to. We present the 
full analysis of these element abundance variations in Paper II and 
show that also these indices can be recovered very well. The Fe in- 
dices Fe4383, Fe4531, Fe5270, Fe5406 are well reproduced by the 
new model, while Fe5015/Fe5335 are slightly too weak/strong in 
the model. Offsets are largest for the two reddest Fe indices Fe5709 
and Fe5782 that clearly cannot be regarded well calibrated (see also 
TMB/K). Finally, the indices Ca4455, NaD, TiOi, and Ti0 2 are 
not sensitive to a/Fe ratio variations and are not well calibrated 
(see also Thomas et al.;^ 2003| >. 

To summarise, the set of indices that turns out to be most 
useful for element abundance ratio studies includes the Balmer 
line indices H5a, H7A, and H7F, the metallic indices CNi, CN2, 
Ca4227, G4300, C 2 4668, Mg 1; Mg 2 , Mgb, and the Fe indices 
Fe4383, Fe4531, Fe5270, Fe5335, and Fe5406. We will use these 
indices to derive element abundance ratios for the globular cluster 
data in Paper II. 

The model data and globular cluster index measurements are 
available at www.icg.port.ac.uk/~thomasd. 
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